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We use phenomenological nonlocal Lagrangians, which lead to non trivial forms for the quark 
propagator, to describe the pion. We define a procedure based on previous studies on non local 
lagrangians for the calculation of the pion parton distributions at low Q 2 . The obtained parton 
distributions fulfill all the wishful properties. Using a convolution approach we incorporate the 
composite character of the constituent quarks. We evolve, using the Renormalization Group, the 
calculated parton distributions to the experimental scale and compare favorably with the data and 
draw conclusions. 
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I. INTRODUCTION 



Parton Distributions Functions (PDFs) and Generalized Parton Distributions (GPDs) 0. 12. Q ■ relating PDFs and 
electromagnetic form factors, encode unique information on the non perturbative hadron structure (for a recent review, 
see 4] ) . Any realistic model of hadron structure should be able to calculate them. We shall proceed in here to study 
the parton distribution functions of the pion in a previous developed formalism 

From the experimental point of view the PDFs and GPDs of the pion are difficult to determine. Pions decay into 
muons or photons and therefore the pion distribution functions can not be obtained from direct DIS experiments. 
The pion parton distribution has been inferred from the Drell-Yang process 0, an d direct photon production [{| 
' in pion-nucleon and pion-nucleus collisions. These set of experiments have been analyzed in ref. ^(j obtaining that 
C^) , the fraction of moment of each valence (anti)quark in the pion is about 0.23 for Q 2 = 4 GeV 2 . 

The parton distribution functions of the pion has been a subject of much discussion in the literature. In 0, 
calculations using quenched lattice QCD, for their lowest moments, were per formed. The lack of sea quarks in 
the approximation implies a greater presence of valence quarks. Its PDFs [13l Il5| and GPDs 

have been also 

calculated in the Nambu-Jona Lasinio (NJL) model 0|- I* 1 the chiral limit, its quark valence distribution is as simple 
as q(x) =6 (x) 9 (1 — x) . Once evolution is taken into account, good agreement is reached between the calculated 
PDFs and the experimental results 0]. The pion PDFs and GPDs have also been calculated in a model with the 
sim ples t pseudoscalar coupling between the pion and the constituent quark fields 0], and in the instanton model, 
[l!^,l20j| . In ref. pj| the chiral limit result of the NJL model for PDFs changes such that q (x) goes to zero for x — > 
and x — ► 1. The pion PDFs have also been calculated in a spectral quark model [l8|. 

Also, the pion PDFs have been calculated in a statistical model, without any dynamical assumption, obtaining 
quite reasonable results [2l| . In [23] the GPDs of the pion are calculated in the bag model. A relevant contribution 
to the calculation of the pion GPDs on the light-front has been given by Tiburzi and Miller |23j , and some remarks on 
the use of the light-front for calculating GPDs can be found in [24|. Looking for a more fundamental approach, not 
related to a unique model, the pion PDF has been also calculated in the framework of Dyson-Schwinger Equations 

m 

In all the previous calculations the pion was built of two valence constituent (anti)quarks. However, these constituent 
(anti)quarks are themselves built of elementary partons. The way to connect a description in terms of constituent 
quarks with the description in terms of the elementary partons was described in |27l l3l[ . Hereafter, we will call this 
procedure ACMP. In |27L 1281 129| the procedure was applied to the nucleon case and in [3(| it was applied to the pion 
case. The same description of the constituent quarks in terms of partons is used for the proton and for the pion. 

The aim of this paper is to study the pion PDFs following the ideas developed in ref 0, 161. Our main objective is to 
develop a calculational method for the PDFs preserving all the fundamental symmetries, i.e. momentum conservation 
and gauge symmetry. Another motivation is to understand the simple result obtained for the pion PDFs in the NJL 
model in connection with our more elaborate description. 
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Our study is done in a lagrangian model [f| built with the following properties: (i) manifest chiral invariance 
for the strong part of the interaction; (ii) a momentum dependence for the quark propagator provided by QCD 
based calculations; (iii) the implementation of local gauge invariance provides the electromagnetic interaction. This 
procedure allows us to have simultaneously a coherent description of the dressed quark propagator, the pion state, 
the dressed quark photon vertex and to preserve gauge symmetry. In this sense our aim is not to fit the best possible 
description of the pion PDFs, but to learn about the influence of the momentum dependence of the quark propagator 
on the definition of the PDFs. The differences between our approach and the previous calculation [2j| will become 
obvious. We achieve our objective in a three step process. The first is to construct the PDFs of the pion in terms 
of its constituents quarks, for that we will use the model developed in ref. 5]. For this model the pion form factor 
was calculated obtaining good agreement with data, and the operator form for the mesonic parton distribution was 
found in terms of the quark propagator and the kernel of the BS equation. The second step is to include the partonic 
structure of the constituent quarks. To do so we will use the ACMP procedure assuming for the quarks the same 
structure functions as was obtained in the analysis for the proton. The last step is to evolve (to next to leading order), 
by means of the renormalization group of QCD, the resulting parton distribution up to the value of Q 2 for which data 
exist pel . 

This paper is organized as follows. In section II we discuss the definition of the PDFs in our formalism and give 
provide some technical details of the calculation. In section III we perform an analysis of the PDFs, with models 
defined from fundamental approaches, achieving a good description of the data. In section IV we calculate the moments 
of the parton distribution, using three different procedures, to check the validity of our method for calculating parton 
distributions, and in section V we provide some conclusions. 



II. THE CONSTITUENT PARTON DISTRIBUTION FUNCTIONS 



In our scheme the pion arises as a solution of the two body BS and therefore is a composite system of two 
"constituent" valence (anti)quarks which determine the PDFs. To calculate them the corresponding operator was 
introduced in ref Q (we refer the reader to this paper for details). 

The parton distribution has three contributions: 



The first one, (x), corresponds to the one body term 
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where Tr represent the trace in Dirac, color and flavor indices, T^ 1 {p,p') is the quark-photon vertex which in general 
can have a complicated structure (see [EUSHi!)- Here T M (p, P) is the BS amplitude for the meson, 



r£ (p,P) = -2, J -0rG a ^s (P,P',P) (is(p' + ip) T M (p',P) is(p'- Ipjj 
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with a, (3, 7 and 8 including spinor, color and flavor indices. The normalization condition for the BS amplitude is: 
2iP» = 
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To calculate the PDF we need the elastic vertex r M (p, p) directly related to the quark propagator by using the 
Ward Identity: 



T»(p,p) 



dS- 1 (p) 
dpu 



(5) 



FIG. 1: Diagrammatic expression of the first one body term in Eq.Q 




FIG. 2: Diagrammatic expression of the second one body term and the two body term in Eq.Q 



The pion momentum can be expressed in terms of the light-front vectors as P M 



4^/2 with p- 



P (1, 0, 0, 1) /\/2, = (1, 0, 0,-1) / (V2-P), where P has a particular value for each frame (for instance P = M/y/2 
for the rest frame). This first contribution corresponds to diagram shown in Fig.^ 
The second term is also a one body term: 
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while the third term is a genuine two body term given by: 



d A P r d 4 P ' 



(2nfJ (2tt) 



is[p'--p) r M (p', p)is[p' + -p 



6(x 



(2 P + P] 



d 



dP u 



+ 5{x 



(2p' + P) 



iS[ P +-P\r M (p,P)tS[p-^P 



d 
dpn 

f3S 



dP u 



c*7 

GafSyS (P',P,P) 



(6) 



(7) 



These two last contributions corresponds to the diagram of Fig [21 and they result from the non local character of the 
interactions . In the contribution q^P (x) the bubble integral has been performed using the BS equation, whereas this 

(2) 

is not possible in the contribution (x) . 

In ref. [5| a particular model was developed defined by means of a chirally invariant lagrangian: 



C(x)=ip (x) (ip- mo) if) (x) + g j| (x) J s (x) + (x) J 5 (x) + g p jj (x) J p (x) 



(8) 
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built with the non local currents defined by 



J s (x) = / d 4 y G (y) ip ( x + ^y] if) ( x - \y 



./.-, I.r) = / ,/ ; G„ (//)<.'( x + iT7 5 ^ (x - iy^ (9) 
J p (x)= I d^yGpiy)^ (x +l-y) -i P tp (x - -y 



2"J r \ 2 
^y\ «f7 5 -0 o 

where u d v = u (dv) — (du) v. The first and second currents require the same Go (y), m order to implement chiral 
invariance, while the third current is self-invariant under chiral transformations. The scalar current, Js, generates 
a momentum dependent mass, and the last current, the "momentum" current, J p , is responsible for the momentum 
dependence of the wave function renormalization. The pseudo-scalar current, J5, generates the pion pole. For 
simplicity, we assume that all the G (y) functions are real. 
Let us define 

G(p) = J d 4 ye^G(y) (10) 

with the normalization condition G (p = 0) = 1. For simplicity we consider that our effective theory arises from the 
large N c limit of QCD. This is equivalent to working in the Hartree approximation, i.e. only direct diagrams will be 
considered. 

The model introduces, in a natural way, momentum dependence to the quark mass and to the wave function 
renormalization of the quark propagator: 



S(p) = Z{p)— . (11) 



where the explicit expressions for m (p) and Z (p) are given in the Appendix [X] The scalar current can give a mass 
for the quark even if the Lagrangian contains no explicit mass term, mo = 0, leading to a spontaneous symmetry 
breaking mechanism similar to that in the Nambu-Jona Lasinio model. On the other hand, the momentum current 
gives rise to a momentum dependent wave function normalization. However, even though the latter contributes to 
the mass, it is not able by itself to break chiral symmetry spontaneously. 
The model defined in Eq. © gives the following interaction terms, 

G a ^ s (p',P, P) = 2 g Go [p') G (p) [s ja S S0 + (ir l5 ) ja [if lb ) 50 ] +2g p G p (p') G p (p) (j>) &fj (12) 

and the corresponding BS equation © can be easily solved leading to: 

T M =^(p,P)= i^T l g vqq G {p) (13) 

In Appendix^ we give the explicit expressions fixing the pion mass and the normalization constant g^qq. This model 
realizes the Goldstone theorem and therefore goes to zero when mo vanishes. This is a consequence of the use of 
the same Go (p) function in the scalar current, Js, and pseudoscalar current, J5, in Eq. (jSJ). On the other hand, we 
have no constrain on G p (p) imposed by chiral symmetry. 
Inserting Eas. l(P2*|> and l(T3|) in Eq. J7J we obtain 

<?i 2 > (x) = 0. (14) 

Summarizing, our scheme leads to two contributions: the standard one, associated with the triangle diagram, and 
defined in @ and a new contribution defined in Eq. I©. This latter contribution arises from the non locality of the 
currents in our model, it is also a triangle diagram, from the point of view of QCD, but once written in terms of the 
BS amplitudes it gets a different structure. 

We now proceed to the explicit calculation of the PDF. First of all we can use the 5-function present in Eqs. (0 
and © in order to integrate over p 3 : 



6(x-~-(2p + P)) =V2Ps(V2p(x-^J -p q ~p a ) (ir>! 
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With this integration, all the quantities appearing in Eqs. (J5J) and (J5J) become linear in p°. In order to show explicitly 
that the results do not depend on the choice of P we use the following variable 



(16) 



The last term has been introduced in order to obtain symmetric expressions. In terms of uj we have 

- m n I - - xj -p ± , 



p = U) \ X 



p — —P 'j = —uj (1 — x) + m\x (1 — x) — p\ 
P+--PJ = ujx + to 2 (1 — x) x — p\ . 



(17) 



Next we perform the Wick rotation in terms of the uj variable and thereafter its integration. 

The calculation is performed in Euclidean space. To bring back these results to Minkowsky space is, in general, a 
non trivial matter, as will be discussed shortly. To check the consistency our scheme we will rely on the verification 
of four fundamental properties of the PDFs which arise form conservation laws, i.e. momentum conservation, gauge 
symmetry, isospin symmetry and the fixed number of valence particles in our state. These properties are: (i) the 
parton distribution must be defined between < x < 1; (ii) the first moment of the parton distribution is 1; (iii) the 
second moment of the parton distribution is 1/2; (iv) the parton distribution must be symmetric q (x) = q (1 — x) . 

The way the first property is realized, in those cases where an exact solution is possible [lfij , is associated to the 
Wick rotation. Let us assume that the masses are momentum independent. The denominators present in the definition 
of the parton distribution, Eqs. (0) and I©, can be written as 



to 2 + ie = — (1 — x) 



2 p 1 , + to 2 
uj — m^x H — = ; i- 



(1-x) (1-x) 
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uj + to 2 (1 — x) — 



pj_ + ni z 



.£ 
i — 

x 



(18a) 
(18b) 



and we observe that the integration over uj is different from zero only if < x < 1. We can perform the Wick rotation 
in the region < x < 1, where this rotation is well defined according to the pole positions given by Eqs. I|18|) . For 
the values of x in the regions x < and x > 1 the Wick rotation is not allowed. 

The models used here are too complex to allow for a similar proof of the first property, but we have been able to 
show numerically that the parton distribution vanishes outside the interval. Moreover, the last three properties come 
out of the calculation and represent our consistency check. 

In many models these four properties are not satisfied. In that case authors proceed in the following way: (i) they 
impose that the PDF is only valid for < x < 1; (ii) they force the normalization condition of the PDF by changing 
the normalization of the wave function or BS amplitude; however in this process, the form factor of the pion will loose 
its normalization; (iii) some authors claim that not satisfying the third and fourth properties is consequence of the 
lack of "gluon component" in the meson in those models. 



III. EXPERIMENTAL ANALYSIS OF PDFS. 



In order to get numerical estimates for the PDFs we need to introduce specific models which provide us with 
expressions for Go (p) and G p (p) or alternatively for m{v) and Z (p) . T here are many papers dedicated to the study 
of the quark mass term in the quark propagator [iM . l.'ilL 133 . I.'ifj. I.'i7t 133. l39l 1 1 (1 |41| . We consider two different 
scenarios based on different philosophies but which give quite coherent results. From now on the analysis is carried 
out in Euclidean space a feature which we indicate by the sublabel E in the corresponding momenta. 

The first scenario, called hereafter SI, is based on the work of Dyakonov and Petrov |3(j, which provides us with 
the momentum dependence of the quark mass term coming from a instanton model. They assume Z (pe) = 1 and 
work in the chiral limit (too = 0). Their results are well described by the expression 

/ A 2 \3/2 

m (p E ) = m + a m [ — 2 (19) 

V A m + Pe / 
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Case 


(qq) 1/a (MeV) 


m (MeV) 


(MeV) 


< r 2 > (fm 2 ) 


SI 


-303. 


2.3 


137. 


0.41 


S2 


-285. 


3.0 


139. 


0.36 


Exp. 


-250 300 


1.5 ~ 8 


135. ~ 140. 


0.44 



TABLE I: Results for < qq > 1 / 3 , mo, the corresponding M v and the rms radius squared < r 2 >, for the two scenarios described 
in the main text. 




FIG. 3: PDFs for SI (left) and S2 (right). We have drawn the contribution from Eq. @ (qi 1} (x), dashed line), from Eq. © 
[Qu^ ( x ), dotted line) and the total parton distribution (full line). 

with A m = 0.767 GeV and a m = 0.343 GeV. 

The second scenario, S2, corresponds to an alternative mass function obtained from lattice calculations by Bowman 

et ai. mm, 

A 3 

m(p E ) = m + a m — m , 5 (20) 

A 3 , + (Pe) ' 

with A m = 0.719 GeV and a m — 0.302 GeV. In their lattice analysis the authors also look for the wave function 
renormalization constant. Their values are reasonably reproduced by: 

/ A 2 \ 5/2 

with a z = -0.5 and A 2 = 1.183 GeV. 

In reference 5] different properties of these two scenarios have been studied. In particular, in tabled we show the 

1/3 

values of (qq) ' , pion mass and mean squared radius obtained for these scenarios. The form factor of the pion is also 
analyzed in obtaining a reasonable agreement. 

In Figs. |3]we show the parton distributions obtained for the two defined scenarios. In both cases we have drawn the 
contribution from Eq.© (qu^ (x), dashed line), from Eq. J^J (q^ (x), dotted line) and the total parton distribution 
(full line). We observe that only the total parton distribution is symmetric around the point x = 0.5, as isospin 
predicts. The first moment of the parton distribution is automatically 1, provided that our pion BS amplitude is well 

normalized using Eq. Q. The contribution to the first moment of the new term, q4 (x), is zero as a consequence of 
charge conjugation symmetry Regarding the second moment of the parton distribution, we have for SI that the 

contribution from q^ (x) is 0.47 and the one coming from <ji (x) is 0.03. In the S2 their values are are 0.46 and 0.04 

respectively. We observe that our "new" contribution, qu (x) , restores the correct value, with a contribution at the 
level of 8%. 

We have also looked to the behavior of the PDFs under variations of the parameters present in the expression of 
m (p) and Z (p). We observe that not all values of the parameters are permitted by our consistency requirements. For 
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FIG. 4: Full line corresponds to S2, Eqs. ^) and The other curves corresponds to changes in the parameters of this 

scenario: (i) az divided by 2 and Az multiplied by 2 (dashed line); (ii) a m divided by 2 and A m multiplied by 2 (dotted line); 
(iii) both az and a m divided by 2 and both A m and Az multiplied by 2 (dashed-dotted line). 

instance, if we increase a m or az by a factor 2 we observe that the second moment of the PDF deviates from 0.5 and 
the PDF is not more symmetric around the point x = 0.5. These results appear because the Wick rotation introduces 
spurious contributions. Let us be more precise, our formalism preserves all the symmetries up to the point in which 
we perform the Wick rotation. It is while performing the latter that, if "spurious" poles and/or cuts are close to the 
physical ones, appreciable errors might arise. The position of these singularities depends on the "models" and the 
"parameter choice" in those models. 

Let us now change our parameters in the opposite direction in order to compare our approach with the NJL model. 
We multiply A m and Az by 2, and at the same time divide a m and az by 2 to insure that the value of the physical 
observables, like the quark condensate, do not change too much. This change makes the currents in Eqs. |0 more 
point like, leading to a model whose results, in the limit -> 0, A TO — > oo resemble those of the NJL model. In Fig. 
0]we show the pion PDF in S2, Eqs. (|20|l and (|21|l . with the following changes: (i) Az multiplied by 2 and az divided 
by 2 (dashed line); (ii) A m multiplied by 2 and a m divided by 2 (dotted line); (iii) both Am and Az multiplied by 2 
and both az and a m divided by 2 (dashed-dotted line). We observe the similarity, especially of the last case, to the 
result of the NJL model in the chiral limit, q (x) = 1. 

Some doubts maybe cast in our calculation having to do with the relation between the calculation in Minkowsky 
and Euclidean spaces, as a consequence of the singularity structure discussed in section II. In section IV we resolve 
this problem by calculating the moments of the parton distribution functions in three different ways and proving that 
for the scenarios described above the three methods fully agree. 

By explicit construction, our quarks are constituent quarks. The experimental parton distributions however, unveil 
the structure of the pion in terms of fundamental partons. How to relate the constituent quarks and the fundamental 
partons is an old problem analyzed by Altarelli, Cabibbo, Maiani and Petronzio |27l |. and applied more recently to 
the pion |3(J and to the nucleon [28ll2^| . Let us recall the main features of this development. 

As shown in ref.j^, the constituent quarks are themselves composite objects whose structure is described by a 
set of functions 4> qoq (x) that specify the number of point-like partons of type q which are present in the constituent 
of type (jo, with fraction x of its total momentum. We will hereafter call these functions, generically, the structure 
functions of the constituent quark. The functions describing the nucleon parton distributions are expressed in terms 
of these new functions 4> qaq (x) and of the calculated constituent density distributions (qo — uq, do) as, 



where q labels the various partons, i.e., valence quarks and antiquarks (u v ,d v ), sea quarks (u s ,d s ,s), sea antiquarks 




9o 



(22) 
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(u s , d s , s) and gluons g. 

The different types and functional forms of the structure functions of the constituent quarks are derived from three 
very natural assumptions j^: (i) The point-like partons are QCD degrees of freedom, i.e. quarks, antiquarks and 
gluons; (ii) Regge behavior for x — > and duality ideas; (iii) invariance under charge conjugation and isospin. 

These considerations define in the case of the valence quarks the following structure function 



i\ ft „\A-l 

2) 

T(A) 



W & Q 2 ) = wTvFT^ - 7= ■ ( 23 ) 



For the sea quarks the corresponding structure function becomes, 

-i2\ & rt ^d- 1 



x,Q 2 hit 1 " 1 ) . (24) 



and, in the case of the gluons, it is taken 



'ft*. , x 



b qog (x,Q 2 ) = -(l-xf- 1 . (25) 



We take the same parameters as in ref.[28j for the proton case A = 0.435, B — 0.378, C = 0.05, D = 2.778, and 
G = 0.135, since these functions should in principle be independent of the hadron under scrutiny. 

Once the convolution defined by Altarelli et al. is performed, we proceed to look for the evolution of the parton 
distribution under the Renormalization Group. The evolution is carried to next to leading order. We find the hadronic 
scale, i.e. the scale at which the constituent quark structure is defined, at Qq = 1 GeV 2 . This scale is determined by 
imposing that the fraction of the total momentum carried by the valence quark at Q 2 = 4 GeV 2 is 0.46 as determined 
from experiment |l0j . 

In Fig. 0we show the evolved PDF for S2 (full line). SI, discussed at the beginning of the section, gives similar 
results. The dashed line corresponds to the evolution of the PDF of S2 without the ACMP convolution. We observe 
that the results are reasonable since not a single parameter was fixed in the fit, although too small in the x ~ 1 region. 
The ACMP procedure pushes the PDF to higher values of x and therefore produces a significant corrections in the 
right direction. 

We show in Fig. [5] also the results obtained with S2 for A m and Az multiplied by 2 and az and a m divided by 2 
(dashed-dotted line) . We observe that the experimental results are well reproduced in this case and our curve is very 
well reproduced by the expression 

xq(x)^ ^ + x a (l-xf with a = 0.62, /3 = 1.04 (26) 

r (a) r (i + p) 

which is quite close to the bests fits in . The fact that the PDF in terms of the constituent quarks is larger in 
the regions x ~ and x ~ 1 leads, after applying the ACMP procedure, to a result which is reasonable in the large 
x region. The NJL model, when the ACMP procedure is applied to it, leads to a PDF which too large in the x ~ 1 
region, although quite reasonable given the simplicity of the model and the fact that no single parameter was fitted 
in the process (small dashed line in Fig. |3J). 



IV. MOMENTS OF THE CONSTITUENT PARTON DISTRIBUTION FUNCTIONS. 

In section II we have calculated the parton distribution q (x) . We have assumed that the parton distribution is 
defined between < x < 1 and have used as a consistency test of this assumption that the first moment of the parton 
distribution is 1 and the second moment is 1/2. From Eq. (JSJ and JTJit can be analitically proved that 

Jdxq(x)=F(0) , 
/ dx x q (x) = \F (0) , ' 

when the integration is not restricted to the region < x < 1 1 . The consistency test proposed in section II reflects 
that the significant contribution to the two integrals l|27|l comes from the physical region < x < 1. 



In 5] is was proved that the pion form factor appears correctly normalized, F (0) = 1, in the formalism. 
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FIG. 5: Full line corresponds to the evolved PDF for S2 including the ACMP procedure. The dashed line corresponds to the 
evolution of the PDF of the same scenario without the ACMP convolution. The dashed-dotted and small-dashed lines are for 
the evolved PDF for S2 with az and a m divided by 2 and A m and Az multiplied by 2 and for the NJL model respectively, in 
both cases the ACMP procedure is implemented. 



The numerical value of the higher moments are not fixed by symmetry arguments. Nevertheless we can use them 
as a proof of the consistency of our procedure if we are able to calculate them in a different way. We proceed to do 
so. 

The moments of the parton distribution are defined as: 



= / dx x s 
Jo 



q(x) 



(28) 



A first procedure to calculate the moments is to perform the integration defined in Eq. I|28[) using the parton 
distribution, q (x) , obtained in section II. An alternative method is to perform the integral over x, using the delta 
functions present in Eqs. (0), JBJ) and Q, before performing the integration over the internal momentum, p, in these 
equations. We focus our discussion in the contribution to (x) s coming from J2J, because the generalization to the 
other two contributions is straightforward. 

Once the integral over x is performed, the contribution to (x) s coming from is 



d 4 p 
(2^ 



(2p + P) 



Tr 



;s( ,,-ipj f M (p,P)iS[p+ l -P 



r„ \v + \p,p + \p)^is(p + \p) v M (,, P) 



(29) 



with the momentum integration restricted to those values that guaranty the condition 



0<--(2p + P)<l . (30) 

The difference between the first and the second procedures is just the order of integration of the variables. Any 
significant problem with singularities must appear in comparing the result of these two methods. 

For calculating the integral present in Eq. (|29|l we need a reference frame. The most adequate one is the infinite 
momentum frame, in which we have = (1, 0, 0,-1) /V2P with P going to infinity. The pion momentum is expressed 
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as 



^TfM'^-ii)- (31) 

In this reference frame the dominant terms appearing in the denominators of the propagators are proportional to P, 
[ (p ± P/2) 2 oc ±-^7| (p° — P 3 ) P ], which lead to an integrand vanishing as P~ 2 for large values of P, except in a 

small region of p° around p 3 . In order to evaluate the contribution of this region it is adequate to make the following 
change of variable 



p Q =p 3 + jp'° (32) 

where m is any mass in the problem, for instance ffi — m (0) . 

Using H32[) in eq. (|30[) we obtain that the integration over p 3 must be performed between —P/(2^/2) and P/(2y / 2) 
and similarly for p°. Then, the contribution coming from 12911 is 



si 



= — lim 

P^oo 



-P/(2V2) P 2?r 



Tr 



iS[p 



ip) T M {p,P)iS[p+\p 



dp 3 

P/(2V2) 27r 

1 



fp± 

(2nf 



P 2 



+ -P,p + ip j »/* iS ( P + \p ) T M ( P . P i 



(33) 



It is a this stage that we assume that the Wick rotation on the variable p'° can be performed without problems. The 
most dangerous term, [n ■ (2p + P) /2] s 1 , is under control for any value of s. It is not difficult to show that this 
integral has a well defined value as P — > oo and that this value is independent of fh in Eq. i|32[l . 

We have applied the same procedure to the other two contributions to (x) , coming form Eqs. 10 and J7J and we 
have calculated the moments of the parton distribution in the two scenarios, SI and S2, defined in section III. We 
have also calculated them directly from Eq. I|28|) with the parton distributions obtained in section III. The numerical 
results coincide in both scenarios for all the moments. 

Summarizing, in here we have calculated the moments of the parton distribution functions for the two scenarios 
defined in section III in two different ways. Firstly, from Eq. (|28|) using the calculated parton distribution obtained 
from the method explained in section II. In this first method we do not need to define a particular frame for the 
calculation. In the second method, we have calculated the moments directly in Euclidean space and in the infinite 
momentum frame. This second procedure has been justified from an analysis in Minkowsky space, by performing the 
Wick rotation in Eq. i|33|) and the equivalent equations for the other two contributions coming from Eqs. © and 
{7J). The methods are based on significantly different calculations and therefore, they provide a test of consistency. 
Those models which violate the check of consistency defined at the end of section II, will also give different numerical 
results for the moments of the parton distribution functions in the two methods. 



V. CONCLUSIONS. 



We have applied a procedure for calculating the PDFs in a framework, based on the Dyson-Schwinger equations, 
first introduced in ref. 5]. The motivation behind our procedure is the preservation of all the conservation laws 
associated with the fundamental symmetries (Poincare covariance, gauge invariance, isospin symmetry and number of 
valence constituents) which implies that before evolution the PDF of the pion built with only two valence (anti) quarks 
has first moment equal to 1, second moment equal to 0.5, is symmetric around the point x — 0.5 and is defined in 
< x < 1. Our procedure, determined by the appropriate implementation of local gauge invariance, is defined in 
terms of an operator, which incorporates additional contributions besides the conventional triangle diagram 0]. The 
method of calculation follows the method used in for the calculation of the GPDs in the NJL model and can be 
easily applied to any model with similar ingredients. 

Our approach preserves all the symmetries up to the moment in which the Wick rotation is performed. The 
particular choice for m (p) and Z (p) becomes crucial at this moment. An arbitrary choice of m (p) or Z (p) can 
have "spurious" singularities which might lead to appreciable symmetry breaking effects. Our expressions for these 
functions, obtained from QCD studies, are consistent, within our numerical precision, with all the symmetries. Note 
the parton distribution vanishes at the support endpoints, x = and x = 1, a feature which occurs in the instanton 
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model [l9j, but not in the NJL model. The new terms in the operator give no contribution to the first moment of the 
PDF, but a significant contribution (~ 6 — 8%) to the second moment and restore the symmetry of the PDF around 
the point x = 0.5. We can approach the results of the NJL model by a change in the parameters. 

In order to check furthermore our procedure we have calculated the moments of the parton distribution in two 
different ways. The consistency of these procedures provides a strong test for our formalism to determine the parton 
distributions. 

In order to describe the experimental data we have realized that the ACMP procedure is necessary. Nevertheless, 
even incorporating the ACMP convolutions, the results with the given parameterization are not spectacular, they 
fall very low for large x. On the other hand, the NJL model, also with ACMP convolutions, behaves oppositely, its 
PDFs is too large in that region. Since we can approach the NJL model in our scheme by increasing A m we do so, 
changing simultaneously a m to keep the low energy properties as close as possible to the experimental ones. The 
improvement is immediate although we have done no fine tuning. Nevertheless, we do not make a big statement about 
this improvement, since there might be other contributions to the PDFs arising from improvements to our lagrangian, 
which are lacking in our formulation. 

Fundamental studies of QCD provide results which lead to effective nonlocal interactions which arc difficult to 
connect with the data. The analysis carried out here shows that the formalism developed in ref.[f| is well suited for 
this purpose. The formalism can be easily generalized to any model, is based on the Dyson-Schwinger equations, and 
therefore contains all the non perturbative input associated with the bound state structure of hadrons, and is suited to 
preserve all the wishful symmetries, as long as we are able to control the passage from Euclidean to Minkowski space- 
times. Besides these general conclusions, we have also shown, in the calculation of the PDFs, that the models thus far 
developed to describe the structure of hadrons at low Q 2 (scenarios SI and S2) require from a convolution formalism 
of the ACMP type to approach the data. Moreover, surprisingly enough, without any change in the parameters of 
these models, the results, once ACMP is incorporated, are qualitatively correct. Finally, with little changes in the 
parameterization we have also shown that one can easily adjust both low energy and high energy observables within 
the same model lagrangian. We have not pretended to do so precisely, since our model is far too crude to be able to 
explain the data, but our calculations shows that this enterprise is feasible. 
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APPENDIX A: DESCRIPTION OF THE QUARK PROPAGATOR AND PION AMPLITUDE. 

In ref 5] it is studied the model described by the lagrangian of Eq. JSJ) . We obtain from the Dyson equation the 
momentum dependence of the mass and wave function renormalization of the propagator of Eq. I|ll|) : 

m + a G (p) 
1 - a p G p (p) 

1 - a p G p [p) 

The constants ao and a p are directly related to couplings constants go and g p : 



a = i8 N c N f90 I ^G (p) 2 Z(p \™<*>) (A3) 
(27r) P - m (P) + *e 



a p = i 8 N c N f9p f -^jG p (p) 2 ^f,^. (A4) 
J (2tt) p 2 -m 2 (p) + te 

The BS equation is solved leading to: 

r M =^ (p, P) = 1 75 t 1 g„ qq G (p) (A5) 



Defining the pseudo-scalar polarization 



n ps (P*) =- i4JV c JV) 



(2nf 



G 2 {p) 



Z( P +IP) Z{p-\P) [\P 2 -p 2 +m{ P + l 1 P)m{p~\P) 



P+\P) - m 2 (p + ±P) + ^e (p - |P) - m 2 (p - \P) + 2 e 







the pion mass is obtained from the relation: 



i = 2. 9o n ps (p 2 







and the normalization constant g-n qq is given by 



1 



9irqq 



an. 



<9P 2 
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